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UNITED STATES DEPARTMENT OF COfS/IMERCE 
Nat;ional Doeanic and Atmospheric Administration 

NATlOrMAL MARIME FISHERIES SERVICE 
Stiver Spring . MD 20310 



FEB":5 2002 



Roy W, Lowe 
Project Leader 
Oregon Coastal Refuges 
U.S. Fish and Wildlife Service 
2127 SEOSU Drive 
Newport, Oregon 97365-5258 



Dear Mr. Lowe: 

Thank you for your letters of August 16, 2001 and January 4, 2002 requesting the "legal 
definition of the term 'accidental mortality' as used by the National Marine Fisheries Service 
(NMFS)," and specifically in Permit No, 782-1532 and the associated Biological Opinion. We 
also appreciate your patience in awaiting our response on this matter. "Accidental mortality" is 
not defined in the Marine Mammal Protection Act (MMPA), Endangered Species Act (ESA), or 
their implementing regulations, hi recognition of the fact that there is an inherent risk of injury 
and mortality associated with performing certain procedures on wild pinnipeds. Permit No. 782- 
1532 allows up to five research-related accidental mortalities per year. As stated in Condition 
B,5 of the Permit, this includes pups that are orphaned and must be euthanized when the mother 
dies as a result of research- The Holder of Permit No. 782-1532 has reported accidental deaths of 
Steller sea hons related to capture and restraint, as required under Condition D (Reporting 
Conditions) of the Permit. 

Consistent with the broad definition offtake" under the MMPA, ESA, and the implementing 
regulations of these acts, the accidental mortality condition is interpreted by the Permits Division 
to include any mortality resulting fi"om the actions or presence of the researchers while 
conducting permit-authorized activities. This includes, but is not limited to: deaths of pinniped 
pups by starvation following abandonment resulting fi*om disturbance to a rookery or the 
research-related death of a lactating female; deaths of marine mammals due to adverse reactions 
to anesthetics or other chemical agents; deaths of marine mammals caused or precipitated by 
infections resulting ft*om intrusive research procedures; deaths of animals due to capture 
myopathy resulting from the stress of capture and handling; and deaths of animals due to serious 
injuries sustained in attempts to escape or evade capture or in response to stampedes or 
aggressive social interactions caused by research activities. 

^. 
Thus, research-related deaths of Steller sea lions under the above circumstances would not be 
unauthorized takes, provided the number of accidental mortahties does not exceed that allowed 
by the Permit. We are currently preparing an Environmental Assessment (EA) on the effects of \ 
research on Steller sea lions. We will include the toxt of the above paragraph as our 
interpretation of accidental mortality in the draft EA and would appreciate your comments. The 
NMFS appreciates your concerns regarding this issue and your efforts to monitor compliance of .^q.^iq^^., , 
research conducted within the National Wildlife Refuge with the provisions of the MMPA and/^ 

^ Primed on Recycled Paper CJV. NO. 05-1 392 (ESH) 

Defendants' Exhibit 2 




^'^IGWOf' 



Case 1 :05-cv-01 392-ESH Document 1 4-5 Filed 01/27/2006 Page 3 of 1 3 

ESA. We are in the process of reviewing the report regarding activities conducted on the Refuge 
under Permit No. 782-1532, and observed pup mortalities on the Refuge, provided by the Holder 
of Permit No. 782-1532, the National Marine Mammal Laboratory, We would appreciate any 
additional information you can provide that might assist us in our review, as well as in 
completion of the EA. If you have any additional questions, please contact Tammy Adams or 
Ruth Johnson, 

Sincerely, 





D. Terbush 
Chief, Permits, Conservation 

and Education Division 
Office of Protected Resources 



cc: Director, NMML, Seattle, WA 

Chief, Endangered Species Division, NMFS Office of Protected Resources 
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ARIATION IN MICROSATELLITES AND mtDNA ACROSS 
THE RANGE OF THE STELLER SEA LION, 
EUMETOPIAS JUBATUS 

Robert G. Trujillo, Thomas R. Loughun, Neil J. Gemmell, John C, Patton, and John W, Bickham-^ 

Department of Wildlife and Fisheries Sciences, Texas A&M University, College Station, TX 77843-2258, USA 

(RGT^JCP^JWB) 

Alaska Fisheries Science Cetiter, National Marine Mammal Laboratory, National Marine Fisheries Service, 

7600 Sand Point Way NE, Seattle, WA 98115, USA (TRL) 

Department of Zoology, University of Canterbury, Private Bag 4800, Christchurch, New Zealand (NJG) 



Genetic variation at 6 nuclear microsatellite loci with bi parental inheritance and the maternal iy inherited 
mitochondrial DNA (mtDNA) was studied at 3 geographic scales (rookeries, regions, and stocks) in Steller sea 
lions {Enmetopias jubatus). Genetic variation was high in both nuclear and mtDNA markers as revealed by a near 
range-wide survey of 21 rookeries. However, population structure was not well defined, and there was no 
obvious phyjogeographic pattern to the distribution of microsatellite alleles. This contrasts with a clear 
phylogeographic pattern revealed by control-region sequences of mtDNA in which 2 well-differentiated stocks, 
eastern and western, are defined as well as 2 distinct groups, Asian and central, in the western stock. Effective 
migration estimates are consistently higher for the nuclear loci than for mtDNA- The difference in patterns 
between the biparen tally and maternally inherited genetic markers can be explained by relatively high male 
dispersal rates and female philopatry, or else there has been insufficient time since populations have been isolated 
for the nuclear loci to have diverged. It is recommended that the presently accepted stock structure be retained for 
management purposes and that further studies be carried out to test the male dispersal hypothesis. 

Key words: geographic vajiation, microsatellites, mitochondria! DNA, Steller sea lion 



The Steller sea lion {Eumetopias juhatus) is distributed 
throughout the North Pacific Ocean and Bering Sea. Rookeries 
occur along the North Pacific Rim from the Kuril Islands; north 
to the Commander Islands; cast through the Aleutian Island 
chain to mainland Alaska; and south to California (Loughlin et 
al. 1987; Fig. 1), This species has experienced a marked decline 
from an estimated 240,000-300,000 individuals in the 1960s 
(Kenyon and Rice 1961) to an estimated 1 16,000 individuals in 
1989 (Loughlin et ah 1992). Population numbers in the United 
States have declined by about 75% over the past 20 years 
(Ccilkins et al. 1999) with most of the decline occurring in the 
western stock. Several studies have investigated the decline of 
Steller sea hons, but the exact cause remains elusive (Loughlin 
1998), Suggested causes include reduced prey availability 
resulting from commercial fisheries (especially walleye pol- 
lock, Theragra chalco gramma., the dominant prey of Steller sea 
lioas), historical sea lion pup harvest, diseases, predation, 
drastic and short-term environmental stochasticity, and envi- 
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ronmental pollution, Pascual and Adkison (1994) concluded 
that a long-term change in the environment or a novel 
catastrophe is the likely cause of the decline. Menick et al. 

(1997) showed a strong positive correlation between a reduction 
in the diversity of the diet and population decline. Calkins et ai. 

(1998) concliided that nutritional stress is the most likely cause 
of the decline. These studies and others indicate the reasons for 
the decline of Steller sea lions are multifaceted. 

Genetic studies of Steller sea lions have been conducted to 
investigate the phylogeographic structure of the populations, to 
identify management units, and to document the potential loss 
of genetic diversity resulting from the decline. Bickham et al, 
(1996, 1998b) found high levels of genetic variation in the 
mitochondrial DNA (mtDNA) of Steller ,sea lions throughout 
their range. They concluded there are 2 genetically distinct 
stocks. The western stock occurs from Prince William Sound 
west to the Kuril Islands and the Sea of Okhotsk, and the 
eastern stock is distributed from southeastern Alaska through 
central California. The western stock was listed as endangered 
under the United States Endangered Species Act in 1997 and is 
declining, while the eastern stock remains listed as threatened 
and has been increasing. It was further suggested that the 
western stock is comprised of an Asian group including ail 
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Asian rookeries and a central group including rookeries in 
the Aleutian Islands, the Gulf of Alaska^ the Bering Sea, and 
Prince Williann Sound (Bickham et al. 1998b)- Another study 
(Bickham et al, 1998a) compared genetic diversity in samples 
;x>m the Gulf of Alaska taken in the 1970s and 1980s to 
samples taken in the 1990s. They concluded that no 
appreciable genetic diversity had been lost despite the -^75% 
decline in numbers. 

Previous genetic studies on pinniped species using micro- 
satellite loci have shown varying amounts of nuclear gene 
variation, measured as the number of alleles present per locus 
and the number of polymoiphic loci present per species. 
Coltman et al. (1996) examined 7 microsatellite loci isolated 
from the harbor seal {Phoca vltulina) in 7 pinniped species {P, 
vitulina^ P. groenlandica, Halle hoer us grypus^ Cystophora 
cnstata^ Arctocephalus gazella^ A, trap kails ^ and A. forsten). 
They found that the number of alleles per locus ranged from 1 
to 7 in these species, Gemmel et al. (1997) examined 20 
microsatellite loci in 18 pinniped species including the Steller 
sea lion and 4 other otariid species. They found 14 poly- 
morphic loci in Steller sea lions. The number of polymorphic 
loci per pinniped species ranged from 3 to 19 with a mean of 
13. 

Allen et al, (1995) illustrated the usefulness of microsatelHte 
markers for investigating population differentiation and 
structuring in gray seals (/f . grypus). They used 8 microsatelHte 
loci to show significant differences in allele frequencies be- 
tween 2 breeding colonies off the northwest coast of Scotland. 
Goodman ( 1 997) studied 7 microsatellite loci from harbor seals 

;m 12 geographic areas in Europe. He found little local 
abstructure but noted that genetic differentiation increased 



with geographic distance. Regional patterns of differentiation at 
the microsatellite ioci were similar to the patterns observed for 
mtDNA (Stanley et al, 1996). Microsatellite markers thus 
appear to be useful in population genetic studies of pinniped 
species, especially in conjunction with mtDNA, but no range- 
wide survey of any broadly distributed species has yet been 
conducted. 

The objectives of this study are to investigate genetic 
stiiicture across the distribution of the Steller sea lion using 6 
microsatellite loci and to compare the phylogeographic 
structure estimated from microsatellite loci to the pattern 
observed with mtDNA. Except for the study by Lidicker et al. 
(1981) on ailozymes, this is the ist study of geographic 
variation using nuclear markers in this species. Nuclear 
(autosomal) loci are inherited biparentally and provide 
a different perspective from the maternally inherited mtDNA. 
The combination of mtDNA and nuclear loci might offer 
a clearer picture of population stmcture and effective migration 
among populations as well as establish the relative contribu- 
tions of males and females to any observed patterns of genetic 
diversity. 



Materials and Methods 

Study area and sample size. — Tissue ,sampie.s were taken as skin 
punches of the hind flipper, preserved in saturated NaCI with 20% 
dimethyl sulfoxide (DMSO), and stored at room temperature (Amos 
and Hoelzel 1991). Samples were taken from pups at their natal 
rookeries ranging from lony Island in the 5ea of Okhotsk east to Si- 
George Reef in northern California (Fig. !). Coordinates and sample 
si7.es are given in Appendix I. Total sample si/e for this study was 
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208 pups from 2i rookeries, representing i2 regions and 3 stocks and/ 
or groups. 

Laboratory mefliod.s.—Toinl genomic DNA was extracted by 
established protocols (Manialis ct al, 1982). Six polymorphic nnlcro- 
^ateiiite loci were amplified by the polymerase chain reactiort (PCR). 
mplification reactions were 25 ).il in volume and consisted of the 
following: OJ-0.5 pg DNA, 2.5 \x\ lOx Taq buffer (0.1 M Tris-HCI, 
pH 8.5, 0.025 M MgCl2, 0.5 M KCI). Z5 |il dNTP mix (2 mM dATP, 
dTTP, dCTP, dGTP, in 0.1 M Tris-HCI, pH 7.9), I ]il of a 10-pM 
solution of each primer, 1 U Taq DNA polymerase, and to volume 
(25 {il) with deionized water. Sequences and fluorescence labels for 
the 6 primer pairs are presented in Appendix II. 

Microsatellite loci were amplilied on a Robocycler thermal cycler 
(Stratagene, La JoHa, California) iis follows: 32 cycles of 95°C for 45 s 
of denaturing, 50°C for 30 s of iuinealing, and TO'^C for 2.5 min of 
extension. The amplifications were visualized on a 1% agarose gel to 
confirm amplification of loci. Amplilication products were analyzed 
using an ABi 377 automated DNA sequencer (Applied Biosystems, 
Foster City, California) using the program Genescan Analysis v. 2.1 
(Applied Biosystems). The program Genotyper v. 1.1 (Applied 
Biosystems) was used to assign allele sizes to the peaks. 

For mitochondrial DNA analysis, a —450 base-pair segment of the 
control region was amplified, and a 238 base-pair segment of this 
amplification was sequenced as described in Bickham et al. (1996). 
Haplotypes were identified based on nucleotide differences and 
assigned a unique letter code. 

Standard genetic variability statistics (number of alleles per locus, 
expected heterozygosity, observed heterozygosity, and allele frequen- 
cies) were calculated using Arlequin v. 2.0 software (Schneider et al. 
2000). Values were calculated for each locus and each population. 
Genotypic phase disequilibrium was estimated to test for locus 
•'^dependence prior to further analysis using GENEPOP v. 3.1 
I lymond and Rousset 1995). Deviation from Hardy- Weinberg 
-fequiiibrium was tested at a = 0.05 using a Markov chain method and 
Fisher's exact test in GENEPOP v. 3,1, To correct for multiple 
simultaneous comparisons, sequential Bonfeixoni corrections were 
applied using a global significance level of a = 0.05 (Rice 1989). 

Analyses were employed to investigate differentiation at 3 geo- 
graphic scales: rookery, regional, and stock or group levels (Fig. 1). 
Paii^ise comparisons at all 3 geographic scales were made for each 
genetic distance measure, Slatkin's linearized Est (Slatkln 1995), 
which is based on an infinite alleles model, was calculated between all 
population pairs using Arlequin v. 2,0. 

Mitochondria! DNA haplotypc frequencies were calculated for all 
groups at the 3 geographical scales described previously, Haplotype 
frequencies were then used to compute Slatkin's linearized Fst for all 
population pairs, and neighbor-joining trees were constructed using 
PAUP V. 4.0b4a (Swofford 2000) for comparison to those generated 
using nuclear microsatellite data. Migration estimates (Nm) were 
calculated from Slatkin's linearized Fst at the 3 geographic scales. 
Significant values for Fst {P < 0.05) estimated from haplotype 
(mtDNA) and allele (microsatellite loci) frequencies were calculated 
by peimutation testing. 
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Table 1,-— Genetic diversity at six microsateMite loci over all 21 
Steljer s^^iion population^ exximined In this sttid^f. Cgiunsri pleadings 
are as follows: n ^ total number of allele;s; As = riieaii liumbef of 
Jll^es jper p^agulation; Mi — overall e^^ected hetero^ygxa^ty^ Mg — 
mmn e^p^ted heterozygo?!^ per pcfjutajtion; il| ^ niean obsejpvred 
hfM^g(3^gity perpc>piiim^^ = overall alfelib size Variance; :% = 

mearr stlJeliie ste valance per popula^jQn^ 
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Results 

Variation at nuclear microsatellite loci among Steller sea 
lion rookeries. — Microsatellite alleles were scored from 6 loci: 
HG63, HG8.10, PVl 1, PVC29, MI lA, and M2B, Each locus 
was polymorphic with the number of alleles ranging from 4 
(HG8.10) to 14 (M2B) per locus. The mean number of alleles 
per rookery ranged from 2.67 (HG6J) to 7,05 (M2B). Overall 
expected heterozygosities ranged from 0.25 (HG6.3) to 0.87 
(M2B). The mean observed heterozygosity per rookery ranged 
from 0.25 (HG6.3) to 0.81 (Ml 1 A), and the overall allelic size 
variance ranged from 1.82 (HG8.K)) to 46,45 (PVl 1; Table !). 

Each lookery was found to be in Hardy- Weinberg equi- 
librium {P > 0.05) at each locus. Tests of genotypic phase dis- 
equilibrium for each locus pair across all lookeries were not 
signilicant {P > 0.354), and thus the genotypes from one locus 
can be assumed to be independent from the genotypes at another 
locus, 

Slatkin's linearized Fst values calculated between all 
pairwise comparisons of rookeries varied from to 0.206. 
Low but significant levels of subdivision (P < 0.05) were 
detected for 75 of 210 paii^ise comparisons (36%), Fst values 
for region pairs ranged from to 0.187 with 23 of 66 (35%) 
significant {P < 0.05). Values at the stock level were 0.002 
between the Asian group and the central group of the western 
stock and 0.007 (significant, P < 0,05, between the eastern stock 
and the central group) between the eastern stock and both the 
central and the Asian groups of the western stock. Neighbor- 
joining trees were constructed from Slatkin's linearized Fst 
values at the rookery level (Fig, 2A) and the regional level 
(Fig. 3A), 

Mitochondrial DNA variation in Steller sea lions. — ^A total of 
55 haplotypes were ob.served among the 194 individuals 
examined. Gene diversity (Fig. 4) ranged from 0.69 (lony 
Island) to I (Atkins and Marmot Islands), from 0.69 (Sea of 



Fic. 2. — Neighbor-joining trees ba.sed on Slatkin's linearized Fst showing relationships among Steller sea lion rookeries based on A) 
microsatellite and B) mtDNA data. Regional designations are as follows: OKH = Sea of Okhotsk; KUR = Kuril Islands; RUS = Ru.ssia; WAl = 
- -stern Aleutian Islands; CAI = central Aleutian Islands; EA! ^ eastern Aleutian Islands; BER = Bering Sea; WGA = western Gulf of Alaska; 
\ = central Gulf of Alaska; PWS = Prince William Sound; SEA = southeiistern Alaska; NCA — northern California. 
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Okhotsk) to 0.95 (southeastern Alaska), and from 0.8 (Asian 
group) to 0.9 (central group). Nucleotide diversity (Fig. 4) 
ranged from 0iK)4 (Lovushki Island) to 0.016 (White Sisters 
Island), from 0.005 (Sea of Okhotsk) to 0.012 (southeastern 
\Iaska and Prince William Sound), and from 0.007 (Asian 
I group) to 0.01 1 (eastem stock). 

Slatkin's linearized Fst was calculated for each possible 
pairwise comparison at the rookery, regional, and stock or 
group levels. Values ranged from to 0.358 among rookeries 
and from to 0.259 among regions. At the stock or group level, 
values were 0.053 between the eastern stock and the Asian 
group, 0,096 between the eastern stock and the centra! group, 
and 0.092 between the central group and the Asian group. At the 
rookeiy level, 33% of the values were significant {P < 0.05), 
and at the regional level 61% of the values were significant. All 
3 estimates were significant at the stock or group level. 

The neighbor-joining tree constructed from Siatkin\s linear- 
ized Fst values for rookeries (Fig. 2B) revealed 3 of the 5 Asian 
group rookeries clustering together at the base of the tree with 
a 4th (Medney Island) branching off immediately. The four 
eastern stock rookeries cluster next, while the central group 
rookeries along with Antsifcrova Island (Asian group) form the 
last major cluster. At the regional level (Fig. 3B), a geographic 
pattern of population subdivision also is clear with clustering of 
the eastern stock regions, clustering of the Asian regions 
(western stock, Asian group), and clustering of the central 
regions (western stock, central group). 

In order to estimate levels of dispersal in Steller sea lions and 
to compare nuclear and mitochondrial gene flow, migration 
\ )timates (Nm) were calculated irom Slatkin's linearized Fst at 
Jthc 3 geographic scales (rookeries, regions, and stocks or 
groups). For nuclear genes, Nm represents the number of 
effective migrants per generation, whereas tor mtDNA, Nm 
represents the effective number of female migrants per 
generation. Values of Nm are greater for nuclear loci than for 
mtDNA in 141 out of 209 comparisons among rookeries and in 
55 out of 65 comparisons among regions. At the stock or group 
level, values of Nm based on mtDNA were 5.41 migrants 
between the central and Asian groups, 5.21 migrants between 
the central group and the eastern stock, and 9.43 migrants 
between the eastern stock and the Asian group. At the stock or 
group level, values of Nm based on nuclear loci were 249.76 
migrants between the central and Asian groups, 74.84 migrants 
between the central group and the eastern stock, and 67.03 
migrants between the ea.stem stock and the Asian group. 



Discussion 

Variation at 6 nuclear microsatcllite loci in Steller sea lions 

was high in term.s of the total number of alleles per locus and 

obsen/ed heterozygosity, but population differentiation was 

'^w. Of the 21 rookeries sampled, 4 (lony Island, Chirikof 

land. Hazy Island, and Ugamak Island) had .signilicantly 
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diflerent values of Slatkin's linearized Fst at the majority of 
pairwise comparisons. At the regional level, 3 had significantly 
different values at the majority of comparisons (Sea of 
Okhotsk, Bering Sea, and northern California). Based on 
mitochondrial DNA, levels of genetic variation were also high, 
but patterns of population subdivision were more defined. At 
the rookery level, significant values of Slatkin's iinearized Fst, 
which indicate significant levels of population differentiation, 
were observed in 36% of the rookery pairwise comparisons for 
microsatcllite loci and 33% of the rookery pairwise compar- 
isons for mtDNA. At the regional level, the percentage of 
significant Fst values remains constant for microsatellites (35% 
of regional pairwise comparisons) and increases substantially 
for mtDNA (61 % of regional pairwise comparisons). The trend 
continues at the stock or group level, where 33% of the Fst 
values are significant, indicating significant levels of popula- 
tion differentiation, for microsateilite loci while 100% of Fst 
values are significant based on mtDNA. 

A distinctive pattern of macrogeographic variation is ob- 
served in the mitochondrial genome, but it is not apparent in 
the nuclear genome. There are at least 3 potential explana- 
tions for this pattern. First, gene flow in Steller sea lions might 
be relatively high but gender biased. This could be explained 
by reference to the breeding behavior and life history of the 
Steller sea lion and the inheritance patterns between mtDNA 
and the nuclear loci. In Steller sea lions, females arc highly 
philopatric to their natal rookery, and males are more likely to 
disperse in search of a breeding territory. Thus, a system of 
male-biased gene flow could explain the lack of geographic 
structure in the nuclear loci of Steller sea lions. In a recently 
published study on dispersal and rookery fidelity by Raum- 
Suryan et al. (2002), no interchange of breeding Steller sea 
lions was observed between the western and eastern stocks 
over a 24-year time period and branding of over 8,500 Steller 
sea lion pups. However, the authors report that an 1 1 -year-old 
adult male that was branded in 1976 on Marmot Island 
(western stock) was observed in the eastern stock during the 
annual breeding season in 1987. They note that the animal 
would have been reproductively and physically mature enough 
to maintain a territory and breed (Raum-Suryan et al. 2002). In 
the mitochondrial genome, maternal philopatry is evident in the 
degree of subdivision s&en at the regional level. Migration 
estimates based on microsatcllite and mtDNA data support the 
conclusion of relatively high female philopati-y and high levels 
of male dispersal. Migration estimates based on microsatellites 
are on average greater than those estimated from mtDNA. 

A system of male-biased gene flow has been hypothesized to 
explain differences in patterns of subdivision in other species. 
High male dispersal with subsequent male-biased gene flow 
has been proposed for DalFs porpoise, Fhocoenoides dalli 
(Escorza-Trevino and Dizon 2000); African elephants, Lox- 
odonta africana (Nyakaana and Arctander 1999); and northern 
populations of the yellow warbler, Dendroica petechia (Gibbs 
et al. 2000). Furthennore, patterns of high variation but low 
differentiauon of microsateilite loci also have been found in the 
African buffalo, Syncerus caffer (Simonsen et al. 1998), and in 
the pied fiycatcher, Ficedula hypoleuca (Haavie et ah 2000). In 
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showing felatibnsSps amoag Stelkr sea lion regions. Tree A is based 
on itiiGrosatellte, and tree B is Mscd on mtDNA data. Abbreviations 
as for Fig- 2. 

species with philopatric females and highly dispersing males, 

the mitochondrial genome is expected to show greater 

-copulation subdivision than is the nuclear genome. Female 

iilopatry helps to preserve any historical patterns of sub- 



division, while male dispersal tends to mask the effects oFsuch 
occurrences in the nuclear genome. 

The. 2nd explanation for the difference in patterns of 
subdivision betv/een nuclear and mtDNA loci is that the 
species population is in fact currently subdivided and effective 
migration among the populations is low. This obviously 
explains the pattern seen in the mtDNA, but why is the same 
pattern not observed in the nuclear loci? When populations arc 
isolated, they will evolve reciprocal monophyiy in approxi- 
mately 4N generations, according to simulation models and 
population genetic theory (Neigel and Avise 1986). Because 
effective population size for mtDNA is only one- fourth that of 
nuclear loci, population divergence will more quickly be 
observed at this locus than for the nuclear loci {Moritz 1994). 
Therefore, apparent extensive gene flow in the nuclear genome 
could simply be the result of insufficient time having elapsed 
since the populations have been isolated to allow for significant 
lineage sorting to have occuired as a result of genetic drift. 

A 3rd explanation tor observed differences in patterns 
between the nuclear and mtDNA markers is lack of resolution. 
This study is based on 6 nuclear loci, all of which were derived 
from phocids. It is possible that studies employing additional 
loci, which should include loci derived from Stcllcr sea lions, 
will confirm or refute the difference in patterns wc propose. 

Notwithstanding the fact that a stronger phyiogeographic 
signal is apparent for mtDNA than for the nuclear micro- 
satellite loci, some degree of population structure exists for 
the nuclear genome as well. In particular, lony Island, the 
only rookery examined in the Sea of Okhotsk, is well 
differentiated by microsatellite loci but somewhat less so with 
mtDNA. The Sea of Okhotsk is at the extreme western edge 
of the range of S teller sea lions, and it might have ex- 
perienced some degree of isolation, either presently or in the 
past. Additional rookeries need to be examined to better 
assess the level of diversification of this region, hi addition, 
microsatellite loci reveal a closer relationship between the 
Asian and central groups of the western stock of S teller sea 
lions than between either of these and the eastcm stock, 
whereas the neighbor-joining tree for mtDNA shows a closer 
relationship between the Asian group and the eastern stock. It 
seems that there is some important phyiogeographic signal in 
the microsatellite loci. 

An argument can be made for recognizing 2 stocks of S teller 
sea lions (eastern stock and western stock — Bickham et aL 1 996) 
and 2 groups (substocks) within the western stock (Asian group 
and central group— Bickham et aL 1998b) based solely on 
mtDNA divergence. Rookeries from the Sea of Okhotsk, the 
Kuril Islands, and Russia (Medney Island and the rookeries on 
the Kamchatka Peninsula) constitute the Asian group of the 
western stock, and rookeries from the Aleutian Islands, Gulf of 
Alaska, Prince William Sound, and the Bering Sea constitute the 
central group. Rookeries from the eastern stock include those 
from southeastern Alaska through California. However, based 
on nuclear data, we find a lesser degree of population structure. 
These 2 views have somewhat different management implica- 
tions. Should these populations be managed as separate stocks or 
as a single unit? Although high levels of male-mediated gene 
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flow might be characteristic of this species, it is not yet proven, 
and these stocks might be reproductively isolated to a greater 
degree than a consideration of the nuclear gene patterns alone 
would suggest. Moreover, the fact that only the western stock is 
declining would argue that there are other biological differences 
between the two populations. Due to the low rates of female 
migrants from the eastern stock to the western stock based on 
mtDNA migration estimates, recovery of the western stock is 
largely dependent on increases in productivity and survival of 
Steller sea hons within the western stock. Because the sample 
sizes used in this study arc not large. In terms of both numbers of 
loci and numbers of individuals sampled, further studies are 
required to test the hypothesis of high rates of male dispersal in 
Stelier sea lions. It thus seems prudent to continue to accept the 
2-stoek hypothesis for purposes of management decisions. 
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Appendix I 

Localities ami sample sizes. — The 208 Stelier sea lion tissue 
samples examined in this study weie taken from pups at their natal 
rookeries from throughout their range. All individuals were analyzed 
for microsatellite variation, and 194 siunples were analyzed for 
mtDNA (italici^ed sample sizes where different). The exact localities, 
regions, and stocks and groups are as follows. Western stock, Asian 
group: Sea of Okhotsk, lony Island 56°24,3'N, I4.3"23'E (// == 9): 
Kuril Islands, Lovushki Island 48'*32.6'N, I53''5(),4'E (// ^ 9), 
Antsivcrova Island 50"ILrN, 154^58.1'E (/; = 10); Russia, 
Zheleznaya Bay 53M.7'N, 160°L3'E (// -^ 10), Medney Island 
54''4rN, I67°40'E {n = 10). Western stock, central group: western 
Aleutian Islands, Buldir Island 52'^1.6'N, I75''58.4'E {n ^ H), 9): 
central Aleutian Islands, Seguam Island 52°2rN, 172^34.4' W {u - 
10), Ogchul Island 52^'>9.7'N, I68°24.2'W (/? ::- 10); eastern Aleutian 
Islands, Bogoslof Island 53°55.7'N, i68^24.2'W (// ^ 10). Uganiak 
Island 54"I3.5'N, 164"47.5'W {n ^ 10); Bering Sea, Walrus island 
57°irN, I69^56'W (n - 14, /2); western Gulf of Alaska, Pinnacle 
Rocks 54°46'N, i61''45.9'W {n = 10), Atkins Island 55^3.2'N, 
1 59^1 7.4' W (/? ^ 5); central Gulf of Alaska, Chirikof Island 
5.5^46.5'N, I55^39.5^W (/? = 9, /O), Mamiot Island 58"13.6'N, 
t5]M7.8'W {n = 10, 5); Prince William Sound, Fish Island 
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59'^52.9'N, I47^20XVW {n =■ i(), 9^ Seal Rocks W9.Wn, 
I46''503'W (n = 1!, 10). Eastern stock: southeastern Alaska, White 
Sisters Islands 57°38'N, l36''i5'W (/? -=^ iO). Hazy Island 55°52'H, 
i 34^34' W (« :- II), Forester Island 54''50'N, 1 33^31' W (/? ^ 11,6); 
northern California, St. George Reef 4r44.7'N, 124^15'W (/? -- 9). 

Appendix II 

Microsatellite primer sequences. —The primer sequences of the six 
microsatellitc loci employed in this study were identified from cloned 
DNA isolated from 3 pinniped species. In this study, these loci were 
found to yield polymorphic PCR products in Steller sea lions. The 
locus name, foiAvard and reverse primer sequences, fluorescent label 
used as a tag for autOvSequence analysis, source species, and primer 
reference arc as indicated: locus M2B: F: 5'-TCAGTCCACCCACC" 
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TAG, R: 5^CCGACTGCTGGGGTAAA(?, TET ftuorescent label, 
Miroiinga angustirostris (Hoelzel et al. 1999); locus PVC29: F: 5'- 
GGTTAATTGTGTTGTTTACATCT, R: 5'-AACCAGAAGAATA- 
GAATTAGCAT, HEX ftuorescent label, Phoca vitulina (Coltman et 
al 1996); locus HG6-3: F: 5'-CAGGGGACCTGAGTGCTTATG, 
R: 5^GACCCAGCATCAGAACTCAAG, FAM fluorescent label, 
Halichaerus grypus (Alien et al. 1995); locus HG8,10: F: 5N- 
AATTCTGAAGCAGCCCAAG, R: 5'-GAATTCTTTTCTAGCA- 
TAGGTTG, TET fluorescent label, H. grypus (Alien et al. 1995); 
locus Mil A: F: 5'-TGTTTCCCAGTTTTACCA, R: 5'-TACATTCA- 
CAAGGCTCAA, FAM fluorescent label, M. angustirostris (Allen 
et aL 1995); locus PVil: F: 5'-GTGCTGGTGAATTAGCCCATTA- 
TAAG, R: 5'-CAGAGTAAGCACCCAAGGAGCAG, HEX fluores- 
cent label, P. vitulina (Goodman 1997), 
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